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Recently, corrections to Einstein-Hilbert action that become important at small curvature are 
proposed. We discuss the first order and second order approximations to the field equations derived 
by the Palatini variational principle. We work out the first and second order Modified Friedmann 
equations and present the upper redshift bounds when these approximations are valid. We show 
that the second order effects can be neglected on the cosmological predictions involving only the 
Hubble parameter itself, e.g. the various cosmological distances, but the second order effects can 
not be neglected in the predictions involving the derivatives of the Hubble parameter. Furthermore, 
the Modified Friedmann equations fit the SNe la data at an acceptable level. 

PACS numbers: 98.80.Bp, 98.65.Dx, 98.80.Es 

1. Introduction 

That our Universe expansion is currently in an accelerating phase now seems well-established. The 
most direct evidence for this is from the measurement of type la supernova .1] . Other indirect evidences 
such as the observations of CMB by the WMAP satellite |2(, 2dF and SDSS also seem supporting 
this. 

But now the mechanism responsible for this acceleration is unclear. Many authors introduce a 
mysterious cosmic fluid called dark energy to explain this. There are now many possibilities of the 
form of dark energy. The simplest possibility is a cosmological constant arising from vacuum energy 
3]. Other possibilities include a dynamical scalar field called quintessence |J, or an exotic perfect 
fluid called Chaplygin gas 

On the other hand, some authors suggest that maybe there does not exist such mysterious dark 
energy, but the observed cosmic acceleration is a signal of our first real lack of understanding of 
gravitational physics 6]. An example is the braneworld theory of Dvali et al. 0. Recently, some 
authors proposed to add a term in the Einstein-Hilbert action to modify the General Relativity 
(GR) |3,|9|- By varying with respect to the metric, this additional term will give fourth order field 
equations. It was shown in their work that this additional term can give accelerating solutions of the 
field equations without the need of introducing dark energy. The matching with observations is also 
discussed for R n model by Capozziello p"(i|. 

Based on this modified action, Vollick used Palatini variational principle to derive the field 
equations. In the Palatini formalism, instead of varying the action only with respect to the metric, 
one views the metric and connection as independent field variables and vary the action with respect 
to them independently. This would give second order field equations. In the original Einstein-Hilbert 
action, this approach gives the same field equations as the metric variation. But now the field equations 
are different from that of from the metric variation. In ref.j^j Dolgov et al. argued that the fourth 
order field equations following from the metric variation suffer serious instability problem. If this is 
indeed the case, the Palatini approach appears even more appealing, because the second order field 
equations following from Palatini variation are free of this sort of instability (see below for details). 
However, the most convincing motivation to take the Palatini formalism seriously is that the field 
equations following from it fit the SNe la data at an acceptable level, see Sec. 3. 

However, the field equations following from the Palatini formalism are too complicated to deal with 
dir ectly . We have to use perturbation expansion to get more amenable approximated equations. In 
ref.[ll|, this has been done up to first order in the perturbation expansion. As will be shown in this 
paper, perturbation approach is valid when redshift is smaller than about 1. Thus when close to 1, 
how good are the first order equations as an approximation to the full field equations? The answer 
we get is rather interesting: the second order effects can be neglected on the cosmological predictions 
involved only the Hubble parameter itself, e.g. the various cosmological distances, but the second 
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order effect can not be neglected in the predictions involved the derivatives of the Hubble parameter, 
e.g. the deceleration parameter, the effective equation of state of the Modified Friedmann equation, 
etc. Furthermore, the first order equations fit the SNe la data at an acceptable level. Thus, we have 
good reason to believe that the full field equations also fit the SNe la data at an acceptable level. 
2. The Field Equations 

First, we briefly review the derivation of the field equations using Palatini variational principle and 
the derivation of the first order approximation to the field equations. For details, see ref.|ll|. 
The field equations follow from the variation in Palatini approach of the action 

S = / d^xyf^LiR) + / d 4 Xy^L M (1) 



2k _ 

where k = 8ttG, L is a function of the scalar curvature R and Lm is the Lagrangian density for 
matter.. 

Varying with respect to g^ v gives 

L\R)Rp V - \L{R)g^ = -kT„ v (2) 
where X^„ is the energy-momentum tensor given by 

In the Palatini formalism, the connection is not associated with g^ u , but with = L'(R)g^ v , 
which is known from varying the action with respect to . Thus the Christoffel symbol with respect 
to h^v is given by 
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r*, = {£„}„ + w [2S^d v) L' - g, v g Xa d a L'] (4) 

where the subscript g signifies that this is the Christoffel symbol with respect to the metric g^ v . 
The Ricci curvature tensor and Ricci scalar is given by 

= R pv (g) - ^{L')- 2 V ^L'V „L' + {L')- X V^ V L' + ^(L'r^V.V^' (5) 



R = R(g) + 3(L') _1 V M V^i' - ^(L')- 2 V ^L'VL' (6) 

where Rfi V {g) is the Ricci tensor with respect to g^ y and R = g^R^v Note by contracting eq.(0), we 
can solve R as a function of T: 

L'(R)R-2L(R) = -kT (7) 

Thus (J5J, iJBJl do define the Ricci tensor with respect to h^. 

Now apply the above Palatini formalism to the special L(R) suggested in ref. 

L(R)=R-^ (8) 

where a is a positive constant with the same dimensions as R and following the factor of 3 is 
introduced to simplify the field equations. 
The field equations follow from eq.© 

a 2 1 a 2 

(i + j^) R ^ - 2 9 ^ R ~m ) = ~ kT ^ (9) 

Contracting the indices gives 



R= l a^L±2Jl + h!£)*\ (10) 
2 a V 4 a 
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Since T is negative and for large \T\ we expect the above to reduce to R = kT, we just take the minus 
sign in the following discussions. 

From eq.@ and ea. (|10|l we can see that the field equations reduce to the Einstein equations if 
\kT\ ^> a. On the other hand, when a ^> \kT\, deviations from the Einstein's theory will be large. 
This is exactly the case we are interested in. We hope it can explain today's cosmic acceleration. 

Recently, Dolgov [12J argued that for a given T MI/ , the field equation for Ricci scalar derived by 
metric variation suffers serious instability problem. Now in the case of Palatini formalism, it can be 
seen from equations J7J) and © that, for a given T^ v , we can directly get R(g) without any need to 
solve differential equations, see ea. (|10|l . Thus there is no instabilities of this sort. This makes the 
Palatini formalism more appealing and worth further investigations. 

Now consider the Robertson- Walker metric describing the cosmological evolution, 

ds 2 = -dt 2 + a{tf(dx 2 + dy 2 + dz 2 ) (11) 

We consider only the spatially flat case, i.e. k = 0, which is now favored by CMB observations and 
also is the prediction of inflation theory . 

We assume a 3> \kT\. In ref.|ll|. it was shown that, up to first order in nT/a, the scalar curvature 
can be obtained from ea. Hl()(l 



Then from eq.©, 



kT 

Rca(-1 + —) (12) 



L^-\a{l-— ) (13) 
3 a 

L'~t(l + —) (14) 

3 V 4a' V ' 

The first order field equation follows from eq.© 

1 3 5 

Rpu = —^ag^v - k(-T^ - j^Tg^) (15) 

The non- vanishing components of the Ricci tensor up to first order follow from eq.(0 and ca. (|lll) 

a 3kT 

Roo = 3- + — 16 

a 8a 



act kT a 2 kT 

Rij = - (aa + 2a 2 + + — — )6 ij 17 

4 a 8 a 



Assume the matter in the recent cosmological times contains only dust with T = —po/a 3 , where po 
is the present energy density of dust. Then 

f = 3p-, f = 3p(--4(-) 2 ) (18) 
a a a 

Now the first order equation describing the evolution of Hubble parameter H = a/a in this modified 
gravity theory, i.e. the Modified Friedmann (MF) equation, can be get from equations ljPo|) . (|T7|) . (fTHfl 
and ifTB}: 
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up 



H = ff^r ( 19 ) 

It is interesting to note that the coefficient a of the R~ l term in the modified Einstein-Hilbert action 
just behaves like a cosmological constant in the recent cosmological times. In ref.Q, it has already 
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FIG. 1: The dependence of the critical redshift ze which characterizes the valid region of perturbation expan- 
sion of the field equations on the matter energy density fraction Q m 

been shown that its magnitude should be about 10 _67 (ey) 2 in order to be consistent with today's 
comic acceleration. This is now easy to interpret. 

Now we must check whether the assumption a ^> \nT\ holds and when this assumption breaks 
down. In order to do this, define ze by a — npo(l + ze) 3 - This parameter gives the upper redshift 
beyond which the assumption breaks down. In terms of ze, ea. (|19fl can be rewritten as 



where as usual, fi m = po/p C Q- Po, Hq is the present matter density and Hubble parameter respectively 
and pco = 3Hq/k. In the following discussions and numerical computations, we always take the value 
Hq = 70kmMpc~ 1 sec~ 1 which is now favored by CMB observations 0- 
Setting z = in cci. (|2"U|l gives 



The dependence of ze on fl m is drawn in Fig.l. In particular, fi m = 0.1,0.2,0.3 gives ze — 
2.351, 1.603, 1.226 respectively. We will use these values in the following discussions. Thus, the first 
order Modified Friedmann equation is at most a good approximation to the full field equation up to, 
e.g. in the case of il m = 0.3, about redshift 1. We want to know up to what redshift first order 
approximation is good enough. This is not so obvious, because when z *~ 1, nT/a ~ 0.73. This is 
not very small, thus the effects of higher order terms in perturbation expansion may not be negligible 
when apply the MF equation on these redshift close to ze- Future high redshift supernova observations 
would reach such high redshift region so the behavior of the MF equation in these region is important. 
On the other hand, if we find that the first order equation is doing good in describing our Universe's 
evolution up to redshift ze and the second order correction is small, this is a strong indication that 




±l(l + zf + ±(l + z E ) 3 
2+!(l^) 3 



(20) 
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the full field equations also behave well in describing the cosmological evolution. In the below we will 
see that the effects of the second correction is rather interesting: the first order equation fits the SNe 
la data in an acceptable level and the second order correction is negligible up to the critical redshift 
ze- But the second order effects can not be neglected in describing the acceleration rate, i.e. the 
deceleration parameter. Thus we have good reason to believe the full field equation is also good in 
fitting the supernova data and this modified gravity theory may be a good candidate for explanation 
of the cosmic acceleration. 

Thus in order to investigate the valid region of the first order approximation, let us consider the 
second order approximation to the theory. 

First, the scalar curvature can be obtained from ea. (|l(Jfl 

+ ---(-)-] (22) 

From eq.JSJ): 

L~a[-- + - -( — ) 2 ] (23) 

3 3a 12 a 

. 4 r kT 1 ,kT^ . 
L -3 [1+ 4a- + 8(-^ < M > 

The second order field equation can be obtained from eq.© 

R nu = -^oig^u - ftf-T^ - ^Tg^] + ~ ^ T 9^A (25) 

The non- vanishing component of Ricci curvature tensor up to second order is obtained from cq. © 

_ 'a 3 kT 15, kT., 9 kT kT . . 

i?oo = 3- + -— + -( — ) 2 + - (26) 

a 8 a 64 a 32 a a 



„ . o aa kT a 2 kT 9a 2 , «X\ 9 3aa kT kT 3a 2 kT kT ^ „ 

R l]= -[aa + 2a 2 + — — + -— + —(— 2 + — + — }6 ij 27 

4a 8a 64 a f6 aa 32 aa 

The second order Modified Friedmann equation now follows from equations (|26|l . (|27|) and l|25|) : 

H 2 = 8«P+32 a + 2 (2g) 
u T 4 a 

The only correction to the first order equation is the quadratic term in the numerator. It is interesting 
to note that while the first order MF equation is formally similar to the Friedmann equation in the 
ACDM model, the second order MF equation is formally similar to the modified Friedmann equation 
in the RS II brane world cosmology model with an effective cosmological constant [l4| . 
In terms of the ze defined above, ea. lffify can be rewritten as 

H = H n m 3, 1+M3 (29) 

Set z — in ea.l(2"5|l gives 



Ho —9^ (Ho — 9^2 _l f3 _ llo "iO 
g *t"m y v 8 m / ' \2 16 m ) m 1 

ZE=( ^/O_llo r ' _1 (30) 

/ 16 m 

The dependence of ze on f2 m is almost identical to the first order equation. The difference is very 
small and is of about order 10 -4 . 
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FIG. 2: The dependence of luminosity distance on redshift computed from the first order modified Friedmann 
equation. The dotted, dashed and solid lines correspond to Q, n — 0.1,0.2,0.3, respectively. The little crosses 
are the observed data 



3. Data Fitting with SNe la Observations 

It is the observations of the SNe la that first reveal our Universe is in an accelerating phase. It is 
still the most important evidence for acceleration and the best discriminator between different models 
to explain the acceleration. Thus, any model attempting to explain the acceleration should fit the 
SNe la data as the basic requirement. 

The supernova observation is essentially a determination of redshift-luminosity distance relationship. 
The luminosity distance is obtained by a standard procedure: 

D t = (l + z)r = (l + 2 )f ^ (31) 

where H(z) is given by ea. H2U|) or ea. l|29|l for first order and second order approximation respectively. 
In the data fitting, we actually compute the quantity, 

M o = 51og(-^) + 25 (32) 
Mpc 

The quality of the fitting is characterized by the parameter: 

* =E g 2 . + g 2 . (33) 



where (Iq A is the observed value, /i o i is the value calculated through the model described above, a 
is the measurement error, af llz i is dispersion in the distance modulus due to the dispersion in galaxy 
redshift caused by peculiar velocities. This quantity will be taken as 

C-mz = q- &z (34) 
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FIG. 3: The difference between the luminosity distances computed from the first order and second order mod- 
ified Friedmann equations. The dotted, dashed and solid lines correspond to Q m = 0.1,0.2,0.3, respectively. 
The difference is very small even up to high redshift larger than 1 which is the upper limit for the validity 
of perturbation expansion of the full field equation. It shows that we can trust the prediction of first order 
equation in calculating luminosity distance in almost all the redshift region smaller than the critical value ze- 

where following ref. Q, a z = 200km/s. We use data listed in ref. [13) which contains 25 SNe la 
observations. Since our purpose is to show that first order approximation fits the data at an acceptable 
level and second order effects can be ignored even up to high redshift such as 1, also because trustable 
observations around redshift 1 is very few, we think 25 low redshift samples is enough. Also we do 
not perform a detailed \ 2 analysis, this is suitable when we get more high redshift samples. 

Fig. 2 shows the prediction of the first order MF equation for fi m = 0.1,0.2,0.3 respectively. 

Fig. 3 draws the difference between the /j,q computed using first and second order MF equations for 
il m = 0.1, 0.2, 0.3 respectively. The prediction of second order MF equation is almost indistinguishable 
from the first order equation up to redshift 1. The difference is about the order 10 -2 even up to redshift 
1. So we did not draw the corresponding luminosity distance curves for second order approximation. 
And this confirms the assertion we made above. We do can trust that the predictions made by the 
first order equation in calculating the luminosity distances are good approximation to the full field 
equations even up to high redshift, just not exceeding the critical redshift ze- 

Fig. 4 shows the dependence of the x 2 computed using ea. l|55|) on the Q m . We can see that \ 2 
gets a little smaller for smaller value of fi m . We think this is an interesting feature of this modified 
gravity theory: It has already provide the possibility of eliminating the necessity for dark energy for 
explanation of cosmic acceleration, now it can be consistent with the SNe la observations without 
the assumption of dark matter. Thus we boldly suggest that maybe this modified gravity theory can 
provide the possibility of eliminating dark matter (There have been some efforts in this direction, see 
ref. [ig). This is surely an interesting thing and worth some investigations. 

However, we should also note that we only use SNe la data smaller than redshift 1 to draw the 
above conclusions. It can easily been seen in Fig. 2 that the differences between predictions drawn 
from the different values of f2 m become larger when redshift is around 1. Thus, future high redshift 
supernova observations may give a more conclusive discrimination between the parameters. 

4. The Deceleration Parameter 




FIG. 5: The dependence of the deceleration parameter q(z) on redshift z for fi m = 0.1,0.2,0.3 respectively. 
The second order effects is seen to be large. 
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Given the observation that our Universe is currently expanding in an accelerating phase, the decel- 
eration parameter should become negative in recent cosmological times. The deceleration parameter 
is defined by q = —d/(aH 2 ). 

In the case of first order approximation, from ea. (|19|l we can get 



21 Kp _ 99(KP\2 
16 a 32 V a > 



H (P. -L. 9 H£. HE -I- l\ 
\ u ~ 4 a A 8 a ^2' 

In the case of second order approximation, from ea. (|28() we can get 



(35) 



_ 16 a J_ 32 V a > J_ 256 \ a > (oc\ 

q ffi-L 9K£WllK£ , 11/K£\2 I I 1 ) 

\ u ' 4 a A 8 a ^ 32 V a > ^2' 

Since up/a = (1 + z) 3 /(l + ze) 3 , we can get the dependence of q{z) on redshift, this is shown in 
Fig. 5. It can be seen that the second order effects can not be neglected for all the three values of f2 m . 

Combined with the result of Sec. 3, we can draw the conclusion that the second order effects can be 
neglected on the cosmological predictions involving only the Hubble parameter itself, e.g. the various 
cosmological distances, but the second order effect can not be neglected in the predictions involving 
the derivatives of the Hubble parameter. This later assertion can also be confirmed by deriving the 
effective equation of state of this modified gravity theory. 

5. The Effective Equation of State 

In this modified gravity theory, following the general framework of Linder and Jenkins |17| . we 
can define the effective Equation of State (EOS) of dark energy. Written in this form, it also has 
the advantage that many parametrization works that have been done for EOS of dark energy can be 
compared with the modified gravity. 

Now following Linder and Jenkins, the additional term in Modified Friedmann equation really 
just describes our ignorance concerning the physical mechanism leading to the observed effect of 
acceleration. Let us take a empirical approach, we just write the Friedmann equation formerly as 

H 2 /H 2 = fl m (l + zf + 5H 2 /H 2 (37) 

where we now encapsulate any modification to the standard Friedmann equation in the last term, 
regardless of its nature. 

Define the effective EOS oj e ff(z) as 

1 d\n8H 2 , . 

^ w ^- 1 + 3^(iT^ < 38 > 

Then in terms of w e //> equation eq.J33 can be written as 

H 2 /H 2 = n m (l + zf + (1 - £l m )e 3 S° (39) 

Using the above formalism, we can compute the effective equation of state corresponding to the 
first and second order approximations. 
The first order EOS: 



The second order EOS 



9 ( 1+z \9 _|_ o( 1+z \6 i 13 I 1+z \3 

M>\l+z E ' ~ oy -l + z E > ' S^I+ze' 

l+z E > S^I+ze' + 2A z +4^i- 

iq 3 / 1+z \9 , / 1 + z \6 I I 1+z \3 

,2 _ -i 16 16\l+z E > ^\l+z E ' T U+z E i 

/ 13/ 1 + z \6 _ 5/ 1 + z \3 , IW9 , 3/ 1+z \3\ 
( 32\l+z E > S ^ 1+ze ' + 2A z +4ll+z E / J 



e ff ~ ( 3/ 1+z \fi 5/ 1+z 13 1 1UO 1 3/ 1+z \3\ 14 ' 

)S: 

13 -(- 

2 = 1__ 16\l+z E ' 1 U+z E < ' \l+z E ' / 11 \ 

e ff 8 (~ll{J±z\6 _ §.( J±2-\3 -L l\0 J- §.{ 1+z \3\ 



These two EOSs are drawn in Fig. 6. The difference is quite obvious. Furthermore, they diverge 
before ze- Concretely, the first order EOSs diverge at roughly z — 1.65, 1.05, 0.75 for fi m = 0.1, 0.2, 0.3 
respectively; the second order EOSs diverge at roughly z — 1.8,1.15,0.85 for tl m = 0.1,0.2,0.3 
respectively. In summary, for higher value of fl m , the divergence is more rapid and the first order 



FIG. 6: The dependence of the effective equation of state u> e // on redshift z in the case of fi m = 0.1, 0.2, 0.3 
respectively. Solid and dotted lines correspond to first and second order MF equations respectively. The 
second order effects is seen to be large. 

EOS diverges more rapid than second order EOS for the same fl m , which indicate that second order 
EOS is a better approximation to the full effective EOS. But we do not know whether the EOS 
computed from the full field equations also undergo such a divergence. Maybe this divergence is only 
an indication of the breakdown of perturbation approximation. Note w e ff < 1, which is also a possible 
behavior of EOS, see ref . (18| . 

6. Discussions and Conclusions. 

Applying Palatini formalism to the modified Einstein-Hilbert action is an interesting approach 
that is worth further investigations. In particular, it is free of the sort of instability encountered in 
the metric variation formalism as argued by Dolgov et. al. |12| . And more importantly, we have 
showed that the MF equations fit the SNe la data at an acceptable level. In order to compute other 
cosmological parameters such as the age of the universe, we need the full MF equation, which has 
been investigated by us elsewhere 0|. It is shown there that the age of the universe is compatible 
with the age of the globular agglomerates observed today. 

We have discussed the effects of second order approximation. We showed that the second order 
effects can be neglected on the cosmological predictions involved only the Hubble parameter itself, 
e.g. the various cosmological distances, but the second order effects are large in the predictions 
involved the derivatives of the Hubble parameter. Thus, in the future work dedicated to discuss the 
cosmological effects or local gravity effects of this modified gravity theory, second order MF equation 
(|28|l seems a better starting point. Of course, whether the effects of third order effects or even the 
non-perturbative effects can not be negligible in redshift smaller than ze for some situations deserve 
further investigations. 

In summary, adding a or the like terms to the Einstein-Hilbert action is an interesting idea, 
which may originate from some Sring/M-theory [2(j, and looks like a possible candidate for the expla- 
nation of recent cosmic expansion acceleration fact. We can see such modifications may accommodate 
the update observational data indicating our Universe expansion is accelerating without introducing 
the mysterious so called Dark Energy. However, so far as we know the update experimental testings 
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to General Relativity Gravity theory within the Solar system |2l| have repeatedly confirmed that 
the GR is right; now for larger cosmic scale if we should modify the conventional gravity theory to 
confront cosmological observations, it is vital to both our basic understanding of the Universe and 
development for fundamental physical theories. 
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